A mesoscopic Resonating Valence Bond system on a triple dot 
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We theoretically introduce a mesoscopic pendulum from a triple dot. The pendulum is fastened 
through a singly-occupied dot (spin qubit). Two other strongly capacitively coupled islands form a 
double-dot charge qubit with one electron in excess oscillating between the two low-energy charge 
states (1,0) and (0, 1). The triple dot is placed between two superconducting leads. Under realistic 
conditions, the main proximity effect stems from the injection of resonating singlet (valence) bonds 
on the triple dot. This gives rise to a Josephson current that is charge- and spin-dependent and, as 
a consequence, exhibits a distinct resonance as a function of the superconducting phase difference. 

PACS numbers: 73.63.Kv, 73.23.Hk, 74.50+r 
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By analogy with quantum optics, the production of 
entangled states in condensed matter devices have in- 
herently emerged as a mainstream in nanoelectronicsi 
Entanglement between electrons, besides checking fun- 
damental quantum properties such as non-locality, could 
be exploited for building logical gates and quantum com- 
munication devices. The realization of electron entangled 
states might result from strong interactions in nanoscopic 
systems. Mostly, this offers a room to treat the spin in a 
quantum dot as a qubit 2 with generally a quite long deco- 
herence time 3 ^*. Spin entanglement scenarios have been 
envisioned in such a framework^. Theoretical research on 
the possibility to control and detect the spin of electrons 
through their charges has also blossomed recentlji§i&. 
Let us recall that the direct coupling of two quantum 
dots by a tunnel junction might be used to create entan- 
glement between spins^ and such spin correlations might 
be observed in transport experiments 9 . Another mech- 
anism to induce spin correlations would consist to place 
a double quantum dot away from resonance in a vertical 
configuration between two superconducting (SC) leadsiS. 
We go beyond and (theoretically) explore other entangle- 
ment mechanisms based on triple dot devices and more 
precisely on the prolific proximity between a spin and a 
double-dot charge qubit. We consider the idea to realize 
singlets resonating between equivalent low-energy con- 
figurations on a triple dot by placing the latter between 
two SC electrodes. This can be viewed as a mesoscopic 
Resonating Valence Bond (RVB) systemii resulting in a 
Josephson current with a distinct resonant-like profile as 
a function of the superconducting phase difference. 

The pillar of this mesoscopic RVB system is what 
we refer to as the mesoscopic pendulum. We take two 
strongly capacitively coupled quantum dots (say, dots 1 
and 3 of Fig. 1) that form a charge qubit; only one elec- 
tron in excess is permitted on dots 1 and 3 which thus 
embodies the weight of the pendulum. The pendulum is 
fixed through the dot 2 which is singly-occupied and off 
resonance (spin qubit). Then, we place the quantum pen- 
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FIG. 1: (Color online) The pendulum with resonating singlet 
bonds under consideration: The dots 1 and 3 are strongly ca- 
pacitively coupled and form a two-level system characterized 
by the degenerate orbital states (1, 0) or (0, 1). Dot 2 is singly 
occupied and the direct tunnel coupling between dots 1 and 2 
(or 2 and 3) is negligible. The spin entanglement induced by 
the SC leads between dot 2 and, say, dot 1 (or dot 3), embod- 
ies the rod of the pendulum. The singlets can also resonate 
through the direct tunnel coupling T between dots 1 and 3. 
Consequences in a SQUID geometry are properly analyzed. 



dulum between two SC leads. The rod of the pendulum 
in Fig. 1 represents the emergence of spin entanglement 
between dot 2 and, say, dot 1, induced by the proxim- 
ity from the SC leads. Through the tunnel coupling T' 
between dots 1 and 3 and the closeness of the SC leads, 
hence the singlet bonds will resonate on the triple dot. 

Model. — Assuming that the capacitive coupling be- 
tween dot 1 (dot 3) and dot 2 is negligible (that is quite 
legitimate for the situation of Fig. 1), the general charg- 
ing energy for the dots 1 and 3 takes the standard form 12 
E n , n , = E cl (n g i - ni) 2 + E c3 (n g3 - n 3 ) 2 + E cl3 (n gl - 



ni)(n g 3 — n 3 ), where we have introduced the charging 
energies E cl<3 = e 2 C E3 ,i/[2(C E iC S3 - C 2 3 )], E cl3 = 
e 2 Ci 3 /(CEiC S3 - Cf 3 ). Here, %i, 3 = C gl , 3 V gly3 /e are 
the mean numbers of holes on the gates — coupled to 
dots 1 and 3 through the capacitances C g i j3 — being 
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used to change the numbers of electrons n\ and on 
dots 1 and 3 through the gate voltages V g \$, C13 stands 
for the capacitive coupling between dots 1 and 3, and 
Csi,3 = Cgi,3 + Ci3- Of interest to us is the strong inter- 
dot capacitive limit C13 ~ Cgi,3 with n g \ ~ n g3 ~ 1/2 
such that £01 ~ £10 < (-^00,-^11) < (£20, £02). The 
low energy physics can thus be studied within the re- 
stricted Hilbert space in which only the (rii, rz.3) = (0, 1) 
and (1,0) states are allowed for the dots 1 and 3. The 
manipulation of a single charge in a double dot is now 
well accessible experimentallj*iiiii^ii&. A similar regime 
can be reached with SC dotsi^i^. We assume that the 
dots 1 and 3 form a nonlocal charge qubit (with only one 
electron in excess) and thus we can resort to the project- 
ing operators P01 and Ao to project on the (0, 1) and 
(1,0) states respectively^. Below, we refer to H13 as the 
charging Hamiltonian for dots 1 and 3 related to E nin3 . 

The dot 2 contains one electron in excess at the energy 
e 2 = — e with e > (spin qubit) and is subject to a 
strong on-site Coulomb repulsion U which is typically 
the charging energy on dot 2, leading to the Hamiltonian 

U 2 = £2 ^ d 2a d 2a + f/t^t^T 4| d 2J. ■ (1) 
a 

The SC leads are described by the BCS Hamiltonian 

U sc= £ / 7T £ 4(#)*.w (2) 

+A i (r)*{(r)*J(r)+H.c, 

where ilj is the volume of lead j, x(r) — (—ifiV + 
| A) 2 /(2m) - fi, and A 3 -(r) = Aj-e - ^W i s the pair po- 
tential. For simplicity, we assume identical leads with the 
same chemical potentials A/x = and = Ar = A. Fi- 
nally, the relevant tunneling Hamiltonians take the form 

H T = J2 ^e i *^- ra) *t(r i )d 2<r + H.c. (3) 

3=L,R;<r=T4 

+ ^ T ie ^- ri )vl/t (ri)dff p 10+ H.c. 

<r=T4 

+ ^ r 3 e^ r «- r ^vI/t( rj ,) dCT p 01+ H.c.. 

<r=T,i 

W T ' = |r'|e a e^ (r3 - ri) Q + + H.c. 

We have introduced the Aharonov-Bohm phase <fi(r n — 
r m ) = J r r ™ d\ ■ A and $ = /ic/2e is the SC flux 
quantum. Here, the symbols 1*1, r 2 , and r 3 refer to the 
positions of the dots 1,2, and 3, respectively whereas 
tl and r^, embody the coordinates in the SC leads L 
and R. Note that in the present setup, an electron on 
dot 1 can either tunnel into the SC lead L or still onto 
dot 3 via T 1 = \T'\ exp(zA). Similarly, an electron on 
dot 3 can either tunnel into the SC lead R or on dot 
1. The tunnel coupling between dot 2 and, say, dot 1 
(or dot 3), is negligible. Since there is a single electron 
in excess delocalized between dots 1 and 3 we find it 



convenient to introduce a unique electron annihilation 
operator d a such that S a d a — 1, d a PiQ = d lcr Pi , 

and d^Poi — cfeo-Poi, with d\ a {dza) annihilating cx- 
plicitly an electron on dot 1 (3). Along the lines of 
Refs. flfllEl. we have also exploited the double-dot charge 
qubit notations for Hq-i. The raising operator Q + 
acting exclusively on the state I711713) — ensures that 
each time an electron travels from dot 1 to dot 3 this 
causes a flip (1,0) -> (0,1) that means Q+| 10) = |01), 
Q-|01) = 1 10), Q z = (P01 - Ao)/2, and Q+|01) = 0. 
Bear in mind that, for H.t = 0, the double-dot charge 
qubit is embodied by the ground-state wave function 12 
|*q) = Vf^llO) - VP^e-^- r ^+ lX \01) with P 10 = 
<* q |Ao|*q) = (Ao) and P 01 = (* q |P 01 |*q) satisfying 
(9q\Ht'\^q) = -\ T '\ whcn p oi = Pio = 1/2. Below, 
we consider T\ = T 2 = T3 = T. The total Hamiltonian 
reads H = H + H T + H r > with H = H 13 + H 2 + H SC - 
Charge qubit at resonance. — Let us start with E m = 
E\q = Eq. The low-energy subspace is made of one local- 
ized electron on dot 2 and one electron delocalized on dots 
1 and 3. Now we derive an effective Hamiltonian 7i e *' 
respecting this reduced Hilbert space. The low-energy 
states are well separated by the superconducting gap A 
as well as the Coulomb repulsions U and E 2 q or Eq 2 . We 
assume that dot 2 is small enough such as U — e ^> e 
to hinder a double occupancy on dot 2; by construction 
(£^20,^02) > (^H)-Boo)- We introduce P the projection 
operator on the lowest-energy subspace and K = (1 — P) 
the projection operator on excited states: 

H eff = p U p + p UK 1 KU p ( 4 ) 

E — KTtK 

An expansion in Tlj- can be performed and we identify 
E with Eq — e, i.e., with the ground-state energy as- 
sociated to Hq. We envision the limit where e/A <C 1, 
nevertheless e > to ensure a trapped electron on dot 2, 
and Ego — Eq A such that we can disregard individ- 
ual qubit (and Kondo) contributions. This realm can be 
achieved experimentally by adjusting the gate voltages of 
the dot 2 and of the double-dot charge qubit. Tempera- 
tures are less than e and (Eqq — Eq) (close to zero). 

The main proximity effects with the SC leads contain- 
ing Pio are the nonlocal Andreev tunneling processes of 
Figs. (2a) and (2b) resulting in a nonlocal spin-entangled 
electron state; the second-order contribution gives a con- 
stant. Assuming Eqq -C Eh we omit the superconduct- 
ing cotunneling involving the transfer of the electron on 
dot 2 to dot 3 through the SC electrode R (involving the 
(1,1) charge state) and then the transfer of the electron 
on dot 1 to dot 2 through the SC lead L; however, this 
would only renormalize |T'|. Other (parasitic) SC co- 
tunneling terms can be completely ignored through the 
prerequisite U — e » e (and (E 2Q ,E 02 ) ^> E 00 ). Be- 
low, we build explicitly the effective TL e ^ in terms of the 
phase Lp = 0z,(r L ) - </>r(t r ) - ^ ( dlls + dl2 ) ' A ( n0 ~ 
tice that the integration from to vr, runs via dot 2 or 
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FIG. 2: (Color online) Main nonlocal Andreev mechanisms 
when the electron on the double-dot charge qubit is initially 
on dot 1 and E±o ~ Eoi = Eq . The numbered arrows indicate 
the direction and the order of the electron transfers. Since e <C 
A and Eoo — Eq <C A, we ignore individual contributions 10 
from dot 2 giving a 7r-junction-- and from the charge qubhi— ■ . 
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FIG. 3: (Color online) Ii/Ij versus 8 — 2nf for different values 
of x = \T'\/J. When T' ~ 0, E eff does not depend on <j> R - 
<f>L resulting in Ii = 0. When |T'| ^> J, the sinus-like profile is 
reminiscent of the double dot case—. For \T'\ = J, when the 
interference |(*j|*q>| is fully destructive (|(*j|*q}| -> 0) 
and sin(# — 2nf) changes sign, an abrupt jump in Ii occurs. 



dots 1 and 3): 

H eff = Ho + H T > +j(s 2 -S- (Ao + Pox) (5) 



j e ^-^i-^) I S 2 • S - -j Q + P w 

Je - iv3+#(ri _r 3 ) g _ 1 j Q-p ol _ 



Remember that 0(ri — 



£ J r ri dli 3 • A depicts 



— f 1 

the Aharonov-Bohm phase accumulated to go from dot 
3 to dot 1, <p — (f>(ri — r 3 ) does not depend on A, and 
by construction J > 0. Here, S = dh^f-dp embod- 
ies the spin of the delocalized electron on dots 1 and 
3. An entanglement occurs between the spin on dot 
2 and the spin S of the delocalized qubit. The pro- 
cess of Fig. (2b) enables the singlets to resonate via 

S 2 • S|10) -> S 2 • SQ+Ao|10> = S 2 • S|01). Since Ht> 
does not involve spin degrees of freedom this does not 
affect the spin entanglement. On the contrary, through 
a flip i\o ~^ Pqi this also stimulates the resonance of the 
singlets. The antiferromagnetic exchange coupling reads 



J(5r) = 



2T 2 



Enn — Er 



sm(kp5r) 
kpSr 



exp(-2<5r/0, (6) 



where T = w\T\ 2 Af(0) with Af{0) being the normal state 
density of states per spin of the leads at the Fermi energy 
ep > A and kj? (hp = |kjr|) being the Fermi momen- 
tum. Note that the lowest-order expansion (4) in 7^7- 
is valid in the limit where T <C e + Eqq — Eq. Here, 
£ = Wf/(7tA) with vp the Fermi velocity represents the 
coherence length of the SC leads. Remember that Sr 
denotes the typical distance between two injected spin- 
entangled electrons in a given SC lead in the case where 
one electron stems from dot 2 and the other from, say, 
dot 1. Thus, to have J{5r) non-zero, Sr should not ex- 
ceed the SC coherence length. For conventional s-wave 
superconductors the coherence length £ is on the order 



of micrometers and thus this poses not severe experimen- 
tal restrictions. The suppression of J(Sr) is only poly- 
nomial oc (1/kpSr) 2 . Exploiting parameters for a two- 
dimensional electron gas (InAs or GaAs) coupled to a 
superconductor (Al or Nb) leads to J ~ 0.05 — 0.5K 10 . 
The ground-state energy (the lowest eigenenergy) for a 
singlet state on the triple dot with (S 2 S — t) = —1 then 
can be evaluated in the charge subspace (|10), 1 01) ) : 



E e *f = En — J—e— (J 2 



IT 



2J\T'\ cos(</? — A)) 



l\2 o T\T'\ \U 1 / 2 

(7) 

We emphasize that the mesoscopic RVB system might be 
detectable through the SQUID-ring setup in Fig. 1; note 
that an excited triplet state on the triple dot leading to 
(S 2 S — i) = does not contribute to the supercurrent. 

SQ UID analysis. — The Josephson current through the 
junction reads I\ = (2nc/^o)dE e ^ /d((j)R— c^l)- When a 
Cooper pair accomplishes a roundtrip along the SQUID 
loop, we get the phase sum rule 27r/ = 27r<£>/<!>o = + 
ip — A where $ is the flux threading the SQUID loop, (p.— 
A is the total (gauge-invariant) phase difference accross 
the triple dot, and 9 ^ is the gauge-invariant phase 
difference accross the auxiliary junction (J'). Thus 



h = -I.J- 



\T'\sm( 



2?r/) 



(J 2 + \T'\ 2 - 2J\T'\ cos(0 - 2tt/)) 



1/2 • 



(8) 



with the critical current Ij = 2eJ/h ~ 5 — 50nA. Note, 
in Fig. 1, I s = h + (2eJ'/h) sin6>. Note also that h in 
Eq. (8) shows a resonance which is a distinct feature of 
this mesoscopic RVB setup. Two limiting cases can be 
distinguished. When |T'| ^> J, the dots 1 and 3 form a 
unique effective grain. Here I± = —Ij sin(0 — 27r/) resem- 
bles the supercurrent through a double dot in a vertical 
configurationiS; the extra minus sign results from the in- 
terplay between Ht' and the Andreev process of Fig. 2b. 
When T 1 = 0, we infer that Ji = (up to negligible con- 
tributions from paths involving only dot 3iS). The dots 
1 and 3 are still coupled through the Andreev process of 
Fig. (2b) that is described by the phase ip — cf>(ri — r^). 



4 



Thus, projecting Eq. (5) on the singlet state of the triple 
dot, the charge qubit is embodied by the ground-state 
wave function = y/Pw\10) + v/T^ie i(ip_ * (ri_r3)) |01) 
with Pio = Pox = 1/2 at resonance. This gives E e ^ = 
(^fj\H e f f\$j) = E Q -e- 2 J. The charge qubit seeks to 
react by minimizing the energy independently of 4>r — 4>l- 

Now, let us discuss the highly-resonating situation 
\T\ = J. Owing to |(*j|*q)| 2 = [1 - cos(9 - 2tt/)]/2, 
i w -sin(6> - 27t/)/(2|(*j|* q )|). For half-integer val- 
ues of / if 9 is around zero, we recover a conventional 
sinus-like behavior h = -Ij\T'\sm{6-2nf)/(J+\T'\) ~ 
sin(6»-27r/)/2. On the other hand, for (0-2tt/) ps 
or integer values of /, we observe that |(*j|\Pq)| — > 0. 
Such a destructive "interference" effect between Ht' and 
the Andreev process of Fig. (2b) results in a marked 
jump in the supercurrent (Fig. 3). The denominator in 
Eq. (8) becomes zero whereas sin(0 — 2irf) changes sign. 
Jumps in I\ only emerge in the highly-resonating regime 
|7~'| = J for the singlets and around 9 — 2irf = 2irn (n 
is an integer) where 2j = —Ijsgn{6 — 2irf — 27m). For 9 
close to zero and $ = 0, the supercurrent yields \ = Ij\ 
this stands for a distinct hallmark of this triple dot setup. 

Charge qubit off resonance. — For a finite energy split- 
ting SE — Eqi — Eio > 0, we must add the extra term 
5EQ Z = 5E(P 01 -P w )/2mH eff . When 5E < (|T'|, J) 
we expect that this does not affect much the preceding 
results. On the contrary, for SE ^> (\T'\, J), to minimize 
energy the double-dot charge qubit will satisfy Pio = 1 
and Pqi = 0. Here, 1 01) becomes a virtual (excited) state 
and this will obviously hinder the resonance of the sin- 
glets on the triple dot. For instance, PHt'P = and 

H eff = H + 5EQ Z +j(s 2 -S- ^ Ao (9) 

- J^Jcos(<^- A) (S 2 -S- J Q-Q+P 10 . 



Since Q~Q + Pw = Pio we check that the spin entan- 
glement now always involves the spin of the dot 1. We 
had to go one step further in the expansion of the effec- 
tive Hamiltonian in Eq. jlj by including a final tunnel 
event from dot 3 to dot 1. There is another allowed pro- 
cess where the electron on dot 1 first tunnels onto dot 3, 
hence the electrons from dots 2 and 3 leave into the SC 
electrode R. Finally, a Cooper pair leaves the SC lead L 
allowing to return to the state |10). The singlet contribu- 
tion to h reads -Ij(2\T'\/SE) sin(0 - 2tt/). By analogy 
to the double-dolsiS, the supercurrent profile is sinus-like. 

Conclusion. — We have theoretically introduced a 
triple dot that is a mounting between a double-dot charge 
qubit and a spin qubit. By placing the triple dot between 
two SC electrodes, a spin entanglement can emerge be- 
tween the spin qubit and the spin of the delocalized elec- 
tron on the double-dot charge qubit. When the charge 
qubit is at resonance, the singlets propagate freely on the 
triple dot owing to the equivalence between the charge 
states (1,0) and (0,1) resulting in a mesoscopic RVB 
system and clear predictions. When the resonance of 
the singlets is maximized, i.e., for |T'| = J, we predict 
jumps in the supercurrent occurring for integer values of 
/ = $/$ if & is around zero. To reach |T'| = J, one 
could tune t he tun nel coupling beween dots 1 and 3 sim- 
ilar to Refs. [Ull Since supercurrents can flow through 
a high- mobility AlGaAs/GaAs two-dimensional electron 
gas^, our proposal might be realized in the future. 
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